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Large  transducer  arrays  are  used  in  active  sonar  systems  to 
emit  powerful,  directional  sound  beams.  The  transducer  head 
velocities  oust  be  known  in  order  to  accurately  determine  the  steady 
state  source  level  and  directivity  patterns  of  these  arrays.  The 
values  for  these  velocities  depend  upon  the  electrical  driving 
conditions,  the  equivalent  circuit  of  the  transducers,  and  the 
geometry  of  the  array.  In  addition,  the  calculation  of  these 
velocities  generally  requires  the  inversion  of  a complex  matrix  of 
order  N,  where  N is  the  number  of  transducers  in  the  array.  For  the 
case  where  the  array  has  no  sysmetry  and  possesses  more  than  56 


O elements,  direct  inversion  of  the  entire  matrix  is  not  practical, 

O and  due  to  storage  limitations  may  not  be  possible.  Several  methods 

} are  presented  here  for  confuting  transducer  velocities  for  very 


large  arrays  when  the  self  and  mutual  Impedance  coefficients  and 
the  driving  forces  are  known  for  all  the  frlfmftnt.fi. 


"PATCH  METHOD"' 
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The  "Patch  Method"  (references  (a)  and  (b))  divides  the  array 
into  overlapping  patches  and  computes  the  velocities  by  doing  the 
matrix  inversion  problem  for  each  patch.  This  method,  which  has 
been  incorporated  into  USL-IBM-704  Program  #0717,  is  based  on  the 
following: 
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Consider  the  Thevenin  equivalent  circuit,  as  depicted  in 
Figure  1,  for  the  jth  element,  in  an  array  of  N transducers. 


Thevenin  Equivalent  Circuit  for  Jth  Element  in  an  Am 


GJ  ia  the  Thevenin  equivalent  driving  force;  Z1  is  the 
Thevenin  equivalent  internal  impedance,  which  will  depend  on  the 
operating  frequency,  tuning  reactor,  source  impedance,  and  com- 
ponents of  the  transducer;  Zrj  is  the  total  radiation  Impedance, 
which  is  discussed  in  reference  (c);  Vj  is  the  transducer  head 
velocity.  According  to  Figure  1 we  can  write, 


If  we  define  Zij  as  the  mutual  radiation  impedance  coefficient 
between  the  ith  and  Jth  elements  and  Zij  as  the  self  radiation 
inpedance  of  the  Jxh  element,  we  can  write, 

N 

OJ  Zij*  Vi,  (2) 

i«T  * Vi 

since  Zrj  Zij  ?j  , (3) 
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A method  for  obtaining  the  Inverse  of  a complex  matrix  in  terns 
of  real  matrices  is  given  in  Appendix  A.  If  we  new  divide  the 
array  into  patches  each  containing  K elements,  we  can  determine 
the  velocities  in  each  succeeding  double  patch  by  using  the  follow- 
ing equations : 

M Mt2K 

OJ  «2ZiJ*  Vi  (known)  + Si  Zlj*  Vi  (unknown)  (7) 

Ui  i*  M+l 


which  can  be  written  as, 

Mt2K 

QJ  = FJ  + ZiJ*  Vi  (unknown) 


Zij*  Vi  (unknown) 


where  M is  the  number  of  velocities  which  have  been  confuted  and 
HJ  is  the  driving  force  minus  the  acoustic  load  from  the  elements 
of  known  velocity. 


Figure  2 

Consider  now  a large  planar  array  as  depicted  in  Figure  2. 

The  array  is  separated  into  patches  each  containing  K elements. 

It  is  necessary  for  the  user  of  this  method  to  have  a computer  and 
program  which  can  invert  a complex  matrix  of  order  2K.  The  procedure 
incorporated  in  Program  #0717  ia  aa  follows: 
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1.  The  coordinates  are  generated  for  each  piston  in  the 
entire  array. 

2.  The  real  and  imaginary  parts  of  the  0 forces  are  generated 
for  all  the  array  elements. 

3a.  The  "patch  coordinates"  XXP(J)  and  YXP(j)  are  established 
by  setting  XKP(J)  - XK(J+M)  and  YKP(J)  - YX(j+M),  for  J-l,  2K. 

3b.  The  "patch  driving  forces"  are  established  by  setting 
GP(J)  » G(J+M). 


3c.  The  FJ  are  computed  for  the  2K  elements  in  the  two  adjacent 
patches  under  consideration.  For  (a+bJ  , M-0;  there  are  no  FJ. 

3d.  The  FJ  is  subtracted  from  GFJ  to  get  KJ. 

3e.  The  Zij*  matrix  and  its  inverse  are  computed  for  the  double 
patch.  After  the  first  double  patch,  there  is  an  option  to  reuse 
the  initial  Zij*  inverse  matrix  or  to  calculate  a new  Zij#matrix 
and  its  inverse  each  time.  The  ZlJ*Wtrix  will  usually  be  the  same 
for  all  the  double  patches  in  the  same  array.  The  values  of  the 
mutual  radiation  impedance  coefficients  depend  upon  the  piston 
separations  and  the  geometry  of  the  array,  and  the  distance  between 
two  elements  in  one  double  patch  will  usually  be  the  same  as  the 
distance  between  two  correspondingly  placed  elements  in  another 
double  patch. 

’(J),  are  confuted  using  the 


3f.  The  "pat 
matrix  equation 

To  compute  the  2X  velocities  in  [a+b]  , the  remaining 
velocities  are  assumed  to  be  aero.  M is  then  set  equal  to  M+K, 
and  the  2K  velocities  in  [B+Cjare  computed  assuming  the  velocities 
in  [A]  are  as  Just  determined  and  the  remaining  velocities  are 
aero.  The  two  sets  of  velocities  on  left  end  column  of  [b]  are 
then  cospared.  If  the  velocities  of  the  elements  tested  do  not 
agree  within  a given  tolerance,  the  program  prints  out  which  pistons 
failed  but  continues  the  calculations  and  testing  for  each  succeed- 
ing double  patch  until  all  of  the  velocities  have  been  computed. 

At  this  point,  if  any  of  the  comparison  tests  failed,  the 
user  of  the  program  has  the  option  to  return  to  step  #3a  and  begin 
recooputing  the  patch  velocities  or  to  continue  with  the  remainder 
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of  the  program.  If  the  choice  la  to  return  to  step  3*1  the 
velocities  In  any  double  patch  are  obtained  by  considering  the 
remaining  velocities  to  be  as  last  coeluted.  This  process  of  re- 
confuting  the  patch  velocities  may  be  repeated  until  all  of  the 
comparison  tests  sore  successfully  passed. 

The  "Batch  Method"  has  been  tested  successfully  far  an  array 
of  378  elements  with  54  elements  in  each  of  13  overlapping 
patches.  The  results  are  given  in  Appendix  B.  Also,  for  mmlTr 
arrays  (N  -54)  the  time  for  the  velocity  computation  could  be 
reduced  by  using  this  method.  For  example,  the  computation  tine 
for  a 54  element  case  was  reduced  from  nine  to  four  minutes  by 
dividing  the  array  into  three  patches  each  containing  18  elements. 

However,  there  are  arrays  for  which  the  "Patch  Method"  in 
its  present  form  fails.  Correct  values  for  the  transducer  hand 
velocities  could  not  be  obtained  for  a $4  element  planar,  unshaded, 
close-packed  array  of  circular  pistons,  with  steering  angles  of 
00  and  600,  Rl  - 0,  X*  - -Xself , and  Ka  - 1.0. 

Here  R1  and  X1  are  the  real  and  imaginary  parts  respectively 
of  Z , Xself  is  the  imaginary  part  of  the  self  radiation  impedance, 
K is  the  wave  number  ( 2 tt/X),  and  a is  the  radius  of  the  pistons. 
It  is  possible  that  the  badly  behaved  cases  may  be  solved  by 
increasing  the  overlap  from  1 patch  to  M patches  (1«M<2)  without 
increasing  the  number  of  elements  in  a double  patch. 


"ITERATIVE  METHOD" 

The  second  method  for  confuting  transducer  head  velocities 
has  been  incorporated  in  USL  IBM-704  Program  #0697  end  does  not 
involve  matrix  inversion.  It  has  a much  siigxLer  computation 
scheme  then  the  "Patch  Method",  The  calculations  can  be  made  by 
a self-correcting  process. 

Consider  once  again  equation  #2. 

H 

Qj  - 5L  ZiJ*  Vi  . 
i-1 

which  can  be  written  as 

N Vi 

GJ  - VJ  IEzi.1*  VJ.  (10) 

i«l 
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Dividing  both  aides  of  equation  #10  by 

_QJ 

VJ-  » „ * Vi 

^ ^ij  VJ  > 
i-1 
or 

£1-  , 

VJnew  * N „ ^ Viold 

^ *3  VJold"  * 

i-1 


N 

^ Zij*  Vi  we  obtain 
i-1  VJ 

(11) 


(12) 


Equation  #12  provides  the  motivation  for  Program  #0697. 


This  procedure  is  as  follows: 

1.  Initial  values  are  chosen  for  the  V's. 

2.  These  values  are  Inserted  into  equation  #12  to  obtain  new 
values  for  the  Vs. 

3.  These  last  values  are  used  to  obtain  still  otter  values 
for  the  V's. 

4.  The  process  may  be  continued  until  it  becomes  stabilised. 

5.  There  is  an  option  to  begin  averaging  over  the  iterations 
at  any  stage  of  the  Iterative  process  or  to  continue  without  averag- 
ing. If  averaging,  the  VJold  are  determined  by  the  following 
equation. 

VJold  « (W)(Vj)+(l-W)(VJavg)>  (13) 

Where  W is  a weighting  function,  the  VJ  are  the  currant  values  for 
the  velocities,  and  the  VJavg  are  the  averages  of  the  velocity 
values  confuted  for  the  J elements.  However,  averaging  may  not 
necessarily  increase  the  rate  of  convergence  or  yield  correct 
solutions  to  cases  which  converge. 


In  cases  where  the  array  permits  one  to  have  confidence  that 
convergence  will  be  sufficiently  rapid,  this  method  may  well  be 
preferred  to  all  others. 


6 
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However,  there  are  array*  for  which  the  iterative  method 
converge*  too  slowly,  converges  to  the  wrong  Unit  or  even  diverges 
Sons  of  the  elements  became  acoustically  blocked  in  the  cases 
where  the  velocity  sequence*  converged  to  the  wrong  limit*  Charts 
Numbers  1-3  show  the  convergence  for  80  element  unsteered  arrays 
with  four  quadrant  sysmtry.  Other  iteration  methods  such  as  the 
on*  introduced  inieferenee  (d),  pp  117-127*  nay  produce  better 
results. 


Reference  (e)  describes  the  iteration  method  incorporated  in 
Program  #0077*  This  program  is  used  for  finding  the  magnitudes  of 
the  V's  and  the  phases  of  the  0's  when  the  phases  of  the  V's  and 
the  magnitudes  of  the  Q'a  are  given.  Convergence  was  much  quicker 
in  Program  #0077*  but  blocked  elements  did  occur  in  rare  instances 
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Chart  Number  1 
Close-  facKed  Afra-s  , 0 


O denote*  convergent  case 

X denotes  divergent  case 

Cases  inside  area  bundaJ  b-  dotted  lines  diverge 
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BORDERING  METHOD1 


The  Bordering  Method  which  at  the  present  time  is  not  incor- 
porated in  any  USL  eouputer  program,  is  described  in  reference  (d), 
pp  105-111.  It  is  basically  a method  for  inverting  the  ZiJ*  matrix 
by  successive  borderings.  Consider  the  matrix  partitioned  In  the 
following  manner : 


which  can  be  written  as 
/ Zn-j  Un  \ 


Where  Zn-j  is  a matrix  of  the  (n-l)th  order  the  inverse  matrix 
of  which  is  considered  to  be  known.  Vn  - (Zni,  .....  Z -> ).  i 


cm  uwpspy 

Huts  r r.s 

CQHiUuI  Konr-25u(34) 


n 


f 


USL  Tech  Memo 
No.  ‘>60-68-66 


where  «-n  • Zrm-Vn  Z’^  UR 


Now,  the  Z inverse  matrix  can  be  obtained  by  constructing  in 
succession 

-I 

zu  zi2  A / zi2  zi3> 


(zu)' 


& *22 


% 'hz  ^3 


9 • • • • 


*31  *32 


33 


where  each  succeeding  one  is  obtained  from  the  proceeding  by 
bordering.  Each  step  of  this  process  is  accomplished  through 
substitution  into  equation  #16. 

The  Bordering  Method  does  not  require  the  storage  of  the  Z1J* 
matrix  although  it  does  require  the  eventual  storage  of  its 
inverse.  Also,  this  method  may  be  used  to  considerably  reduce  the 
time  for  the  velocity  computation  in  cases  where  the  Z^j*  inverse 
matrix  is  known  for  an  array  which  differs  from  the  array  under 
study  only  in  the  respect  that  the  former  has  fewer  elements. 
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APPENDIX  A 


The  inversion  of  a complex  matrix  only  requires  a method  of  in 
verting  a real  matrix.  If  we  are  given  a complex  matrix  £ ■ R + Jx 
its  inverse  can  be  shown  to  be 

[H  + JX]-1  » C + JD 

where 

C - Ql  + XR"1  X]  _1 

and 

D « - CXR-1 

Therefore,  to  invert  a complex  matrix,  one  needs  to  invert  two  real 
matrices . 


APPENDIX  D 


The  following  373  element  array  cook  two  hours  and  five  minutes 
of  IBM-704  Computer  time  to  solve  fur  the  transducer  head  velocities 
Two  passes  were  made.  None  of  the  comparison  tests  hod  failed  on 
the  second  pass. 


